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20.  (Contd) 

Meads  immediately  to  the  absolute  potentials  of  the  rocket  sections.  Approxi¬ 
mate  formulae  are  developed  to  calculate  the  sheath  potential  from  the  charge 
in  the  sheath.  A  fast  converging  reiteration  gives  the  charge  on  each 
conductor,  the  sheath  charge,  and  the  part  of  the  conductor  potentials  that 
are  due  to  the  sheath.  Explanations  are  given  for:  (1)  the  charge  on  the  rear 
section  having  the  opposite  sign  of,  and  a  smaller  magnitude  than  that  on  the 
forward  section;  (2)  the  magnitude  of  the  probe  voltmeter  reading  being  less 
than  that  of  the  rear  section  voltmeter  for  moderate  beam  currents,  and  the 
same  for  beam  currents  way  beyond  saturation;  (3)  the  ratio  of  these  two 
voltmeter  readings  approaching  1  with  increasing  magnitude  of  beam  current; 
(4)  the  potential  of  the  probe  being  negative  for  ion  emission;  (5)  the  magni¬ 
tude  of  the  probe  potential  reaching  a  maximum,  then  approaching  7ero  with 
increasing  magnitude  of  beam  current.  T wo  different  methods  are  given  for 
calculating  the  approximate  potential  ajKthe  probe  with  the  probe  removed. 
Calculations  of  charging  time  constant^  show  them  to  be  considerably  larger 
than  the  plasma  periods.  A  few  observations  are  made  on  the  effect  of 
vehicle  size  and  fraction  of  vehicle  surfac^  covered  by  a  dielectric  on  sheath 
thickness,  vehicle  potential,  and  charge. 


1.  Cohen,  H.  A.  ,  Sherman,  C.  .  and  Mullen,  ( E.  G.  ( 197 9 1  Spacecraft  charging 
due  to  positive  emission:  an  experimental  study,  Geophys.  Res,  Letts. 
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Potentials  and  Charges  on  Conducting  Rocket  Sections 


1.  INTRODUCTION 

Rockets  and  satellites  may  become  charged  to  high  enough  potentials  to  present 
a  significant  nroblem.  A  method  of  controlling  the  potential  is  to  emit  an  ion  or 
electron  beam.  Cohen  et  al1  flew  a  rocket  with  three  conductors:  forward  section  (1), 
rear  section  (2),  and  extended  probe  (3)  (see  Figure  1  and  Table  1)  to  258  km 
altitude  at  night.  Conductors  1  and  2  were  separated  by  an  insulating  ring.  At 
certain  times,  an  ion  or  electron  beam  of  known  current,  Ib>  and  particle  energy, 

E^,  was  emitted  from  the  front  of  the  forward  section  causing  certain  steady  state 
potentials  and  charges  to  develop  on  the  three  conductors.  Voltmeters  measured 
=  Vg  -  Vj  and  0^  -  -  V  j.  The  measurements  in  Table  2  are  part  of  the  results. 

They  are  taken  from  Table  1  and  the  two  extreme  points  of  Figure  5  of  Cohen  et  al.  * 
They  were  measured  after  steady  state  was  reached,  so  except  for  secondary 
electron  current,  the  return  current  equals  1^  and  electrostatic  theory  is  used.  An 
ion  beam  was  emitted  with  the  electron  beam,  but  its  current  was  negligibly  small; 
thus  it  is  ignored  in  the  last  line  of  Table  2.  Questions  arose:  What  were  the 
absolute  values  of  the  potentials?  Why  was  ~  0b?  Why  was  0(  <  0^? 

(Received  for  publication  9  November  1982) 

1.  Cohen,  H.A.,  Sherman,  C.  ,  and  Mullen,  EG.  ( 1979)Spacecraft  charging  due  to 
positive  ion  emission:  an  experimental  study,  Geophys.  Res.  Letts. 

6:515. 
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al  a2  a3  bl  b2 


b  12  S 


0.19  0.19  0.07  1.27  0.205  1.475  0.26  1.7  1 


The  purpose  of  this  report  is  to  make  an  approximate  analysis  of  the  measurements 
in  Table  2  and  to  gain  an  understanding  of  some  of  the  basic  electrical  phenomena, 
including  answers  to  these  questions. 

First,  approximate  values  of  coefficients  of  potential,  capacitance,  and  induct¬ 
ance  are  calculated  in  Section  2.  They  provide  a  more  complete  physical  under¬ 
standing,  including  (1)  the  approximate  potential  the  position  of  the  probe  would  have 
without  the  probe,  and  (2)  the  charging  time  constant.  Then  the  sheath  radii,  abso¬ 
lute  potentials  of  the  conductors,  part  of  the  latter  due  to  the  sheath,  charges  on 
the  conductors,  sheath  charge,  charging  time  constant,  and  ion  traverse  time  through 
the  sheath  are  calculated  and  interpreted  in  Section  3.  The  summary  and  conclu¬ 
sions  are  given  in  Section  4  and  five  auxiliary  detailed  calculations  in  Sections  Al 
to  A5  of  the  Appendix. 
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Table  2.  Measured  Beam  Current,  Energy,  and 
Conductor  Potential  Differences 


Beam 

!b 

piA 

Eb 

keV 

*b 

kV 

♦t 

kV 

ion 

1 

0.  2 

0.  05 

0.  02 

ion 

9 

1.  0 

0.  42 

0.  25 

ion 

11 

2.0 

0.  47 

ion 

12 

2.0 

0.  55 

0.  34 

ion 

374 

2.  0 

1.  04 

electron 

-10, 000 

0.  09 

-0.  09 

-0.  09 

2.  COEFFICIENTS  OF  POTENTIAL,  CAPACITY,  AND  INDUCTION 


Let  p.j,  c...  and  c.^  for  i  i  j  be  the  coefficients  of  potential,  capacity,  and 

induction  respectively.  Let  the  subscripts  1,  2,  and  3  designate  the  forward  section, 

the  rear  section,  and  the  extended  probe  respectively.  From  electrostatics,  for 

2 

example,  a  slight  generalization  of  Page, 


V  =  PQ  +  V  , 
s 


(1) 


Q  =  C(V  -  V  )  , 
s 


(2) 


V1 

rv  i 

v  Is 

V  = 

V2 

,  V  = 
s 

V2s 

.V3. 

_V3s 

"pll 

p12 

p13 

P  = 

p12 

P22 

p23 

_p13 

p23 

P33 

V 

Q  = 

°*2 

. 

C11 

°  12 

c  13 

C  5 

C12 

C22 

c23 

c  13 

c23 

c  33 

(3) 
(4),  (5) 


(6) 


(7) 


(8) 


2.  Page,  L.  (1935)  Introduction  to  Theoretical  Physics,  Seconded.,  %  120. 
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P  and  C  are  symmetric,  the  p./s  and  c..'s  are  positive,  and  the  c./s  for  i  4  j  are 
negative.  V\g  is  the  part  of  V.  due  only  to  the  sheath. 

Definitions 

Cj  Capacitance  of  the  forward  section  by  itself. 

Cj'  Capacitance  of  the  forward  section  with  the  rear 
section  attached,  insulated,  and  uncharged. 

Cjg  Capacitance  of  the  forward  and  rear  sections  attached 
and  electrically  connected, 
a  Radius  of  a  cylinder  or  a  sphere, 
b  Half  length  of  a  cylinder. 

It  can  be  shown  that 

ci<cr<c12.  (9> 

3 

Smythe  gave  the  following  empirical  formula  (here  slightly  modified  to  te 
it  good  to  0.30%  instead  of  0.41%)  for  the  capacitance  of  a  right  circular  cyli  il 
conductor  of  radius  £  and  half  length  b: 


c(a.b)  =  [  0.0707  +  0.  0615  (-)  ]  a  nf 

a 


a  and  b  in  meters.  In  particular,  this  gives 


C.  =  0.  0629.  nf  (11) 

1  4 

C12  =  0.  0689Q  nf  .  (12) 

Let  us  carry  out  three  gedanken  experiments,  designated  with  the  subscripts 
a,  b.  and  c.  Vg  =  0  and  Eqs.  (1),  (4),  (6),  and  (7)  are  used  for  each.  In  the  first 
experiment,  charge  q^  is  put  on  conductor  1  while  2  and  3  are  left  uncharged.  Then 

V la  =  Pllqla  (13) 

V2a  =  p12qla  <H> 

V3a  =  p13qla  '  <15) 


3.  Smythe,  W.R.  ( 1962)' Charged  right  circular  cylinder,  J.  Appl.  Phys.  33:2966. 


10 


‘la 
'  la 


II  v  12  V 


la 


13  V 


3a 

la 


liv  dividing  l\qs.  (14)  anri  (15)  bv  Kq.  (13),  this  becomes 


P12 


lY  =  <:n  +  rl2  pTT  *  ° 


^13 

Pjl  13  Pn 


rotating  tho  p.  's  and  ‘s  to  C  j'.  Multiplving  this  bv  p ^  j  and  using  Eq.  (17)  give 
•he  equation  for  i lie*  top  left  element  of  tilt*  product  CP.  )  Physical  consideration 


shows  that  !  V , 


P12  "  P11  • 


\  ,  ,  so  Kgs.  (Id)  arid  (Id)  show  that 

la 


(18) 


In  tho  second  experiment,  charge  q^  is  put  on  1  and  q,  on  2  such  that 
V  ^  =  ^  2b  w’t*1  no  charge  on  3.  Then 

V lb  pllqlb  p12q2b 
Vlb  =  P12q 1 h  ‘  p22q2b 
Vdh  1  Pl:iqib  +  P2dr|2i,  ' 


(19) 

(20) 

(21) 


-1 


Again.  3  is  small  and  far  enough  aw  av  for  V  ](j  to  equal  (q^  +  cl2b*C'12  *  accurately. 

“  (dlt)  +  q2b'f'P'  ^  an<l  q  lb  ^  /'l<l;b  are  suhs,'tule<J  Hlto  Kqs.  (19)  and  (20).  The 
former  may  then  be  written 


P12  "  +  1)C12  *  A  P1 1  ’ 


(22) 
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Substituting  this  into  the  latter  leads  to 


p22=  /‘2p1i-<a2-i)ci2'1.  <23> 

The  value  of  A  is  obtained  in  Section  A1  and  is  Eq.  (A10), 

A.  =  3.  770  .  (24) 


Thus,  Pjj,  p)2,  and  p22  may  be  obtained  given  C^'. 
F rom  Eq.  (15), 


The  value  *  would  have  if  2  were  missing  is  shown  in  Section  A2  and  is  Eq.  (A  1 5 ), 
P130_1  =  °- 2659  nf  •  (26) 


Since  3  is  small  and  far  enough  away,  if  1  and  2  were  electrically  connected, 
P13  1  would  be 


-1  _ 

PI2,3 


qJb 


+  q 

3b 


2b 


(27) 


=  0.  264j  nf  , 

Eq.  (A  16),  as  also  shown  in  Section  A2.  (Note  the  end  of  Section  A2.  )  It  is  reason¬ 
able  to  assume  that 


-1 


130 


'13 


'12,3 


(28) 


similar  to  Eq.  (9).  (Actually,  P , 30‘  *  <  P  , 3'  *  <  Pj 2,  3’  **  >  V3b  =  (q lb  +  q2b>Pl  2,  3 
and  q^  =  A  q  are  substituted  into  Eq.  (21),  yielding 


p23=  (A+  l)Pi2,3  pj3  .  (29) 

in  the  third  experiment,  charge  q3c  is  put  on  3  with  no  charge  on  1  or  2.  Then 
V3c  =  P33q3c  ‘  (30) 
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A  maximum  effect  calculation  shows  that  the  presence  of  1  and  2  changes  V.  by 
less  than  M'’u,  so  this  effect  is  neglected.  Therefore,  since  the  capacitance  of  3  is 

a.,,  q.,  =  a.,V  g  1  in  R1KS  units,  whee  g  =  8.98755  m/nf,  so  that  from  the  value 

a  a  C  a  a  C 

of  a.^  in  Table  i. 


bid  "  1JH-q 

llius,  all  p.  .'s  mav  be  obtained  given  and  p,„  *. 

'll  1  1  > 

ignoring  M,  a  good  approximation, 

>>22  _  .  P12  ''ll 

C1  1  1)  ‘  1  12  C21  1)  ’  e22  TT  ’ 

U  "  >’llt,22  "  P12~  r  ^  +  irt  12  'lC  l'  '  *  4  1 2  '*  • 


(2  1) 


(M2),  (MM),  (M4) 


<  M  5 1 


Choosing  q.,  =  -  q  ^ .  tne  capacit'  between  I  and  2  is 


V— V 


°1  1  ^  22  '  *•  12 
C  12 


1  1 


->>12 "  P22  '  ‘  12°  ' 


(Mti) 


So,  as  the  gap  between  1  and  2  approaches  0,  this  capacity  and  I  c .  - 1  become  very 
large,  since  (.' j  ’  approaches  Cj.,.  Including  M,  as  C'j'  approaches  (j.,.  Id  approaches 
0  and  p  ^  approaches  p  ( .j. 

With  Eqs.  (9)  and  (28)  in  mind,  the  following  reasonable  pair  of  values  of 
effective  length  of  conductor  1  duo  to  2  being  present  and  of  p  ^  ‘  are  chosen: 


C , 1  =  c(a  ,  b  +b„/2)  =  0.  0859.  nf, 
1  112  :> 


(M7) 


where  c(a,b)  is  given  in  Kq.  (10).  and 


-  1 


1M 


1,  -I 

_  2  P  I  MO 


12,  M 


1  =  0.  28 5 0  nf 


(M8) 


These  result  from  using  the  values  in  Table  1,  Kq.  (26),  and  Kq.  (27).  The  value 
in  Kq.  (M7)  accurately  equals  the  interpolated  value,  Ty(Cj  +  Cj9),  from  Eqs.  (11) 
and  (12).  f  rom  the  above. 


15.  16 
12.  07 

M.  77 


12.  07 
2M.  74 


M.  7  7 
M.  8M 

5 

1  28.  4 


(ill)'1 


( M  9 ) 


1M 


M.  8M_ 


From  Eqs.  (3)  and  (39), 


111.  1 

-56.  2 

-1.59 

-56.  2 

70.7 

-0.  462 

pf  . 

(40) 

-  1.  59 

-  0.  462 

7.  85 

Since  Eqs.  (32)  through  (35)  show  that  the  value  of  can  be  critical,  the 
question  naturally  arises,  "How  reasonable  is  Eq.  (37)?".  Rocket  sections  1 
and  2  were  separated  bv  a  two-in.  cylinder  of  fiberglass  which  was  essentially 
hollow  except  for  a  half-inch  thick  fiberglass  disc.  To  check,  the  capacitance 
between  the  two  ends  of  an  18.  7 -cm  length  of  the  rocket,  which  included  the  gap, 
was  measured  and  found  to  be  44.  7  pf.  A  correction  due  to  the  missing  lengths 
of  1  and  2  was  calculated  to  be  10.  2  pf.  probably  slightly  too  large.  Adding  this 

I 

and  equating  to  Eq.  (36),  U  is  obtained  which  from  Eq.  (35)  leads  to  C ^  -  0.  0653  nf. 
not  far  from  the  value  in  Eq.  (37). 


3.  NUMERICAL  RESULTS 

Sheath  radii  are  calculated  from  Eqs.  (A22)  and  (A23)  in  Section  A3  and  are 
shown  in  row  3  of  Table  3.  Using  1^  =  1^(0.  7)  (see  the  end  of  Section  A3),  larger 
values  of  K._  are  obtained;  1.4!)  m  for  1^  -  1  pA,  3.  5%  larger  than  4.  87  m  for  9  pA. 
and  smaller  percentages  larger  than  the  values  tabulated  for  the  other  beam  currents. 
Note  that  the  probe  is  inside  of  the  sheath  except  for  the  1  pA  case. 

Values  of  conductor  potentials  consistent  with  Table  2  are  given  in  rows  4  to  6 
of  Table  3.  V v  is  set  =  0  for  reasons  given  below.  The  values  of  V lor  the  11  pA 

“  I  " 

cases  are  guesses.  The  37  4  pA  case  is  seen  to  be  near  or  at  saturation,  the 
minimum  value  of  |I^I  at  which  |V  =  E^,  so  V  ^  is  chosen  =  -K().  The  total  net 
charge  in  the  sheath  is 

qs  =  "  <ql  '  q2  4  'V  •  l'41) 

V  ,  Q,  and  qfi  are  obtained  for  each  value  of  1^  bv  successive  approximation.  First, 

V  is  guessed.  This,  V  from  Table  3,  Eqs.  (2)  and  (40)  give  1J,  and  Eq.  (4  1)  gives 
q  V  is  then  calculated  from  Eqs.  (A38)  and  (A3!1)  in  Section  A4.  Eqs.  (A37) 
and  (A32)  are  used  for  V.^  for  ■=  1  pA;  Eqs.  (A35)  and  (A3  II  are  used  for  Hie 
other  cases.  The  value  of  b  chosen  in  these  equations  is  b  -  b^  *  b,  2  =  ].  3725  m. 
These  calculations  are  repeated  until  self-consistent  values  are  obtained.  This 
procedure  converges  fairly  rapidly  for  1()  =  1  p A  and  quite  rapidlv  for  the  other 
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cases  since,  for  them,  Iv.l  <  '  1Y1.  The  results  are  shown  in  rows  7  to  13  of 
Table  3. 


l'able  3.  Sheath  Radii,  Conductor  and  Sheath  Potentials,  Charges,  Probe  Potentials 
and  Currents,  and  Charging  Time  Constants 


Row 

1 

lb 

1 

<) 

1  1 

!  1 

12 

374 

-  10,  000  4 A 

2 

Kb 

0.  2 

1.  0 

2.0 

2.  0 

2.  0 

2.  0 

0.  09  keV 

3 

H 

s 

1.  25 

4.  87 

5.  43 

5.  43 

5.  70 

33.  03 

11.  56  m 

4 

v. 

-50 

-420 

-57  0 

'67  0 

-550 

-2000 

90  V 

5 

^  2 

0 

0 

0 

0 

0 

0 

0 

6 

V3 

-30 

-170 

-  100 

-200 

-210 

-  960 

0 

7 

Vls 

36 

67 

80 

95 

75 

46 

-  6 

8 

V2s 

36 

67 

80 

95 

75 

46 

-  6 

9 

V3s 

20 

65 

77 

92 

73 

46 

-  6 

10 

ql 

-7.4 

-50 

-67 

-79 

-65 

-  223 

10.3  nC 

1  l 

q2 

2.  3 

23 

31 

36 

30 

112 

-5.  0 

12 

^3 

O 

l 

-1.  04 

-0.  32 

-1.03 

-1.  19 

-  4.62 

-0.  11 

13 

5.4 

28 

37 

44 

36 

116 

-5.  2 

14 

V  33 

-30.  7 

-134 

-41 

-132 

-153 

-  593 

-14  V 

Id 

V3  ‘  V33 

0.  7 

-  36 

-59 

-68 

-  57 

-  367 

14 

16 

!t 

0.  67 

8.  3 

16 

16 

1 1.  3 

35 

-  3  nA 

17 

Xr30 

0.  93 

0.  91 

0.  91 

0.  90 

0.  84 

-183 

18 

T 

7.  5 

5.  6 

6.  2 

7.  2 

5.4 

0.60 

0.  0010  ms 

19 

tt<T  ) 

0.  9 

2.  2 

2.  4 

9.4 

0.  5 

20 

vv 

13.  8 

1.42 

The  most  striking  aspect  of  these  results  is  that  the  sign  of  qj  is  always  oppo¬ 
site  to,  and  that  of  q9  the  same  as,  that  of  the  beam  particles.  The  reasons  are 
clear.  The  beam  removes  charge  of  one  sign,  leaving  a  charge  (q^  after  reaching 
steadv  state)  of  the  opposite  sign.  After  the  beam  is  turned  on,  the  potential  of 
and  charge  on  the  forward  section  increase  in  magnitude.  Its  charge  tends  to 
induce  a  potential  of  the  same  sign  on  the  rear  section.  Both  sections  therefore 
attract  and  collect  return  plasma  particles  of  the  same  sign  as  the  beam  particles, 
thus  retarding  the  rise  of  the  magnitude  of  Vj  and  V9.  Since  this  current  to  the 
rear  section  is  the  only  current  to  or  from  it  except  voltmeter,  secondary  electron 
current,  ana,  for  art  ion  beam,  secondary  and  emitted  electrons  from  the  probe, 
which  are  negligible  or  small,  this  section,  2,  gains  a  charge  of  the  same  sign  as 
that  of  the  beam  particles.  This  transient  current  remains  appreciable  until  Vj 
stops  changing,  that  is,  until  steadv  state  is  reached.  It  flows  so  as  to  keep  \’ 
near  zero,  that  is,  so  that  V99  nearly  cancels  \  9  j .  The  potential  on  conductor  i 
due  to  charge  q ^  is  designated  V.^  p.^q ^ .  f  or  positive  ion  (electron)  beam  emission. 

V.,  cannot  become  appreciably  positive  (negative)  because  the  return  particles  lack 
sufficient  energy  to  cause  this.  Since  the  rear  section  2,  is  close  to  1  and  almost 
of  comparable  area,  the  magnitude  of  the  negative  (positive)  potential  of  2  required 
to  collect  enough  return  current  to  balance  the  voltmeter  current  must  be  quite 
small— so  it  is  neglected  here,  that  is,  must  be  close  to  zero;  here,  it  is  chosen 
equal  to  zero.  Thus  the  rear  end  of  the  sheath  surface  must  come  in  close  to  the 
rear  section.  The  potentials  constant  and  return  current  collected  by  2  equal  to 

voltmeter  current,  d,  / Z,  ,  plus  current  to  the  probe  emitter,  if  anv,  and  that  due 
b  b  4 

to  secondaries  from  1,  constitute  steady  state.  The  voltmeter  resistance, 

zb  -  109  U  t  5%.  Also,  I  q^  |  <  I  q  j  I .  Physically,  this  is  due  to  1  and  2  not  being 

close  enough  to  each  other  to  be  100%  coupled  and  to  2  being  smaller  than  1.  Mathe¬ 


matically,  this  is  because  V,,^  nearly  equals  -V9  , 
is  small  and  far  enough  away  for  iv  ^1  and 


and  Pj.,  <  p.,.^.  The  probe,  2, 
to  be  small,  less  than  2%  of  |v’j|. 


Basically,  the  same  reasoning  holds  for  the  extended  probe,  so  ~  How¬ 


ever,  since  its  collecting  area  is  much  smaller  and  it  is  farther  from  the  front 


section  than  the  rear  section  is,  with  an  ion  beam  not  much  beyond  that  for  satura- 


lion  (  £  1  m A ) ,  the  voltmeter  current.  lt  0^/ Z(,  is  appreciable.  The  voltmeter 
resistance,  4  Z  -  3  X  I010  Q  i  5%,  Even  with  the  probe  emitter  and  secondary 
electron  current,  Ip,  to  the  rear  section  or  outside  of  the  sheath,  the  probe's  steadv 
state  potential  must  be  appreciably  negative  to  attract  enough  return  ion  current. 


I  ^  to  equal  I(  -  lp.  Secondaries  from  1  to  3  are  probably  insignificant  and  are 
neglected.  The  potentials  constant  and  1^  =  l^  *  1  constitute  steadv  state.  There¬ 
fore,  <  0^.  If  q.j  -  =  0-  Vg  is  essentially  the  probeless  potential,  that  is. 


4.  Huber,  W.  B.  private  communication. 
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saturation  t  n  it  Oie  I ,  *  •<  potential  a  it!  to  be  rjrae'ir  11  '•  /<■)•<  lor  the  probe  to  '.-ollii  *  i 

no  more  return  electron  current  tnar  required  to  huluiv.  e  -1,  t  1 wnero  1  is  tt.i  ! 

probe  emitter  and  second  trv  eleet  it»n  current  *o  the  larvcirti  section.  'I  bus,  for  £ 

triis  i.,se.  .,  aid  1  were  negligible.  V\  lien  q„  -  V  ,  .  l>.  tie-  probe  neither  attracts  > 

nor  rci>els  j»u '•Moles.  but,  tt(  rt  rough  fi  rs*  anproxim’Mori,  p-.ssihlv  collects  1^ 

times  tpe  ratio  of  the  average  solid  angle  the  nrobe  subVnds  at  an  element  of  the  ^ 

snout!:  surface  to  2.;  steridians,  7 

t; 

'.4  2' 

The  probe  radius,  a  =  0.07  m.  The  displacement.  d  -  1.7  1  m,  of  the  probe  fivci. 
the  center  of  the  sheath  makes  only  4  7  error  in  I.  -.q  lor  1  -  f*  (i A  and  less  error 
for  higher  values  of  i  1.  .  Tor  comparison,  values  ol  l4  and  1^  ^  lrom  lap.  (42) 

ire  given  in  rows  Hi  and  17  of  Table  at.  1  is  nearly  independent  of  1^  since,  1 

from  llq.  (,\20i,  l!“  is  nearly  propoiMonal  *o  I  .  I  =  J4  -  I  so  is  'ess  than  . 

Since  1  is  appreciably  less  than  1,  except  possible  lor  the  electron  beam  cast  , 

I  ■  I  consistent  with  the  ~it»0s‘.  of  the  values  of  \  row  '4,  Taiih  .1. 

A  second,  different  (not  equivalent*  wav  of  calculating  tin  approx;.:.  ,te  protv  - 
less  potential  is  as  follows.  (.'.ilculnte  q  and  V  for  q  ,  -  0.  Substitute  it  and  j 

q  ,  =  0  in  tne  bottom  line  of  l,q.  (2!  mil  solve  for  \  ...  l  or  exampb  ,  for  ij  -  ’  m  i  | 

d  fiA,  this  procedure  gives  0.  06  and  i!>  \  as  compared  ♦(>  0.  7  and  -'16  \  respe.  4r  i  -  j 

1\  in  row  15  ol  Table  “■  . 

If  /.  were  increased  from  zero  to  infinity  md  other  given  quantities  remain'-  . 
constant,  the  (steady  state!  value  o|  q.,  would  go  from  a  maximum  magnitude  .aid 
tne  same  sign  >.s  through  zero  to  a  value  of  sm  diet  magnitude  and  opposite  sign; 

V  (  would  go  from  \  an-1  1  iron-  .  :  avin-uir.  to  pntctu  allv  Zero;  lor  an  ion  lv  an 
!  ,,  a nd  1  would  go  from  a  m a  xi mum  '  o  zero;  lie  an  ole  ■!  run  bean .  ]  veuil d  go  i  ron 
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values  of  |  lbl  approaching  saturation  potentials  and  increasing,  the  (constant) 
value  of  becomes  effectively  more  nearly  infinite  (1^  more  nearly  negligible). 

This  is  due  to  1^  becoming  a  smaller  fraction  of  t lie  total  return  current  with 
increasing  |  lbl  .  This  is  borne  out  by  the  ratio  of  ob  to  in  Table  2  being  closer 
to  1  the  larger  |lb(  is.  For  ^  =  374  pA,  tliis  ratio  from  the  values  in  Table  3  is 
inconsistently  large,  indicating  that  the  values  of  Vj  and  V  shown  there  are  chosen 
somewhat  too  negative,  and  this  case  was  below  saturation.  As  |lbl  increases 
from  zero,  I  V.^l  increases  to  a  maximum  near  saturation,  then  approaches  zero 
as  |  lbl  approaches  infinity.  What  makes  the  maximum  ''  Beyond  saturation,  V 
remains  practically  constant,  the  return  current  increases  because  it  includes 
the  additional  beam  current,  so  |  V.(|  must  decrease  to  increase  I  1  |  .  For  an  ion 
(electron)  beam,  q.^  and  \  must  increase  (decrease),  eventually  clianging  signs, 
and  the  fraction  of  lb  that  the  probe  collects  must  decrease  (decrease)  in  order  to 
balance  tile  probe  currents.  The  values  in  row  6  of  Table  3  are  consistent  with 
this  behavior. 

Now  t he  effect  of  secondary  electrons  is  considered.  Let  6  be  the  average 
number  of  secondary  electrons  per  primary  incident  particle.  The  rocket  surface 

4 

was  iridite  treated  aluminum.  For  ion  beams  below  saturation,  values  of  6  for 

N()+  on  this  surface  for  energies  up  to  2  keV  are  needed.  In  lieu  of  this,  ttie  value 

used  is  6  =  0.  18,  ttie  yield  from  1  keV  0+  ions  on  molybdenum.'  The  return 

current,  ignoring  the  small  (•  fv’o  lb)  voltmeter  and  emitter  currents,  even  if  only 

a  fraction  i>  ('2/3)ofthe  secondaries  from  1  exit  the  sheath,  is  l(tj  -  -t  6);  so 

the  sheath  radii  in  Table  3  are  reduced  8  to  1  l°o.  Also,  the  secondaries  decrease 

the  charge  density  in  the  vicinity  of  the  rocket  thus  decreasing  q  and  increasing 

q j  (decreasing  Iqjl  ).  Vj  may  increase  slightly  except  at  or  above  saturation.  The 

other  1  -  v  of  the  secondaries  goes  to  conductor  2.  About  0.  11^  must  return  to  2 

to  balance  this  '0.051  secondary,  '0.051  2  to  1  voltmeter  current,  and  -0.041 

r  r  r 

to  0  current  due  to  secondaries  from  2.  A  more  complicated,  complete  theory  is 
required  for  a  more  detailed  and  accurate  description. 

The  amount  and  effect  of  secondary  electrons  in  the  electron  beam  case  is  much 

6  5 

different.  Using  the  formula  of  Sternglass  and  the  values  for  aluminum: 

6  =  0.  97  and  U  =  300  V,  the  value  6  =  0.72  is  obtained  for  U  =  90  V. 

max  max 

Thus  many  secondary  electrons  must  be  generated  by  the  return  electrons  hitting 
the  forward  section.  Since  they  have  no  more  than  a  few  eV  energy,  how'ever,  they 


5.  St  anna  rd,  1’.  It.  ,  Katz,  I.,  Mandril,  M.  J.  ,  Cassidy,  J.  J.  ,  Barks,  D.  1-J.  , 

Rotenberg,  M.  ,  Steen,  B.G.  (1980)  Analysis  of  the  charging  of  the 
SCATMA  (Bi’8-2)  satellite.  Report  NASA  C R-  1 6534 8,  SSS- K- 8  1 -4  7  98. 
pp  22-24,  29. 

6.  Sternglass,  K.  .1.  (1954)  Rackscattering  of  kilovolt  electrons  from  solids. 

Bins.  Rev.  95:345. 
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do  not  escape,  but,  in  steady  state,  must  be  collected  by  the  forward  section  as 
fast  as  they  are  generated.  Therefore,  the  sheath  radius  is  not  changed  bv  the 
production  of  secondaries,  and  the  return  electron  current  (including  returning 
beam  electrons  if  above  saturation)  equals  the  beam  current.  If  above  saturation, 

Vj  cannot  be  affected.  The  only  essential  effects  of  the  secondaries  are  (1)  to 
increase  the  electron  density  near  the  forward  section,  causing  the  electric  field 
to  be  larger  than  otherwise  near  the  forward  section  and  smaller  in  the  outer  part 
of  the  sheath,  and  (2)  to  increase  q(  and  thus  |  q^.1  .  Other  electrical  quantities 
may  change  slightly.  Some  of  the  return  electrons  will  scatter  instead  of  being 
collected,  causing  the  sheath  radius  to  increase.  This  is  believed  to  be  a  small 
effect  with  this  rocket  and  is  not  treated  here. 

Ionization  by  beam  and  return  ions  is  negligible  due  to  the  small  ion  speeds  and 
high  ionization  threshold  energy.  For  the  1^  =  -10  mA  case,  however,  both  the 
beam  and  return  electrons  produce  ionization  near  the  rocket.  This  has  the  effect 
of  decreasing  q.,  increasing  the  charge  density  (decreasing  | p I )  near  the  rocket, 
and  thus  of  decreasing  lqgi-  Vj  is  unaffected  since  ll^l  is  much  greater  than  the 
saturation  value.  Ionization  at  140  -  150  km  altitude  has  a  moderate  effect  at  most1 
since  the  minimum  ionization  mean  free  path  there  is  about  1000  to  2000  m. 

The  charging  time  constant,  T  ,  of  the  forward  section  is  calculated  and  com- 

3  -3 

pared  to  the  plasma  periods.  Fertile  charge  density,  p.q  =  10  cm  ',  the  latter 

are  0.  76  ms  and  3.5  /us  for  ions  and  electrons  respectively.  The  former  may  be 

calculated  approximately  as  follows.  For  simplicity,  the  probe  is  ignored  since  it 

has  only  a  small  effect,  the  small  voltmeter  current  is  neglected,  and  the  ion  transit 

time  is  assumed  to  be  much  smaller  than  T  .  Let  the  instantaneous  and  final  values 

respectively  of  conductor  minus  sheath  potential  at  conductor  1  be  x  and  X  =  V  j  -  V  j 

at  conductor  2  be  -v„  and  and  of  forward  section  return  current  be  i  ,  and 

2s  2s  rl 

l^j.  From  the  top  line  of  Eq.  (2), 


*1  =  CHX  '  C 


1 2v2s  ‘ 


=  i  ,  +  6i  -  1 


W~  -  lrl  "  rl  ‘b  "  ‘■liar  ‘-12  3T 


"  c  , x  -s  • 


The  assumption  is  made  that 


lrl  X  *rl 


xand  v2s  =  X  V2s 


(45),  (46) 


7.  Leadon,  R.  E.  ,  Woods,  A.  J.  ,  Wenaas,  E.  P. .  and  Klein,  II.  H.  (1981) 
Analytical  investigation  of  Emitting  Probes  in  an  Ionized  Plasma, 


Then 


dv2s  V2s  dx 

3! -  -  3T  nd  |cll  '  c  1 2 


2s  ,  dx 

ar 


=  'lb(1  *  jt*  • 


Assuming  x  =  0  at  t  =  0  when  the  beam  is  turned  on,  the  solution  of  Eq. 
_  t 

x  =  X I 


;[l  -  e'  T  . 


C  (V  “V  )+C  V 
ll'  1  Is'  12  2s 


(47),  (48) 

(48)  is 


(49) 

(50) 


Values  of  T  from  Eq.  (50)  are  given  in  row  18  of  Table  3.  Thus,  the  charging  time 
constant  is  large  compared  to  the  plasma  period  for  ion  beam  currents,  1^  -  12  pA, 
but  not  for  the  largest  ion  beam  current  or  the  electron  beam  current.  As  a  check, 
the  approximate  time,  t^,  for  an  NO4  ion  to  traverse  the  sheath  radially,  starting 
at  t  =  0,  is  calculated  for  each  beam  current  except  11  p A  (see  Section  A5).  For 
ion  beams,  the  ion  is  assumed  to  start  at  the  sheath  boundary  with  velocity  v^  in 
the  -r  direction;  for  the  electron  beam,  the  ion  is  assumed  to  start  at  the  conductor 
with  zero  velocity.  Values  of  t^  for  charging  time  constants  T'  =  t  and  t^  >  T  are 
shown  in  rows  19  and  20  respectively  of  Table  3.  The  charging  time  constant  due 
to  capacitances  is  seen  to  be  much  larger  than  the  ion  transit  time  and  thus  solely 
determines  the  time  for  steady  state  to  be  reached  only  for  the  1  pA  case.  The 
time  to  reach  steady  state  is  determined  by  both  the  capacitances  and  the  ion 
transit  time  for  the  9  and  12  pA  cases,  and  essentially  solely  by  the  ion  transit 
time  for  the  374  pA  and  -10  mA  cases. 


4.  SUMMARY  AND  CONCLUSIONS 

An  analysis  is  made  of  measurements  of  potential  differences  between  pairs 
of  three  conductors  (two  rocket  sections  and  an  extended  probe)  flown  in  the  night¬ 
time  ionosphere  when  an  ion  or  electron  beam  was  emitted  from  the  forward 
section.  1 

Approximate  values  of  coefficients  of  potential  (and  thus  coefficients  of  capa¬ 
city  and  induction)  can  be  calculated  without  using  complicated  potential  theory  for 
the  conductors  of  an  object  which  has  a  fairly  simple  geometry.  These  coefficients 
are  calculated  for  the  forward  section  ( 1),  rear  section  (2)  and  extended  probe  (3)  of 
the  rocket  described  in  Cohen  et  al.  1  Reasonable  choices  of  the  limited  values  of 
two  parameters  are  sufficient  to  give  the  value  of  the  P  matrix  which  is  inverted 


to  give  the  C  matrix.  A  measurement  of  capacity  between  the  forward  and  rear 
sections  confirms  that  the  important  parameter  is  chosen  approximately  correctly. 
The  magnitudes  of  c13,  C23>  and  c^^  are  much  smaller  than  those  of  the  other 
c.  s  since  the  probe  was  small  and  far  enough  away  from  the  other  conductors  to 
interact  with  them  only  weakly.  P  could  also  be  measured  as  follows.  Hang  the 
rocket  away  from  other  conductors  >  10  m  above  ground  and  run  a  small  ground 
wire,  0,  up  to  the  rocket.  Measure  the  capacities  between  1  and  0,  2  and  0,  3  and 
0,  1  and  2,  2  and  3,  and  3  and  1.  The  first  three  measurements  give  the  p^'s, 
nd  the  last  three  give  the  p.Js  for  i  4  j.  The  effect  of  using  more  accurate  values 
of  the  c.  ,'s  than  those  used  in  this  report  is  expected  to  be  small. 

A  formula  is  developed  for  a  cylindrical  sheath  radius  and  for  interpolation 
between  it  and  the  spherical  sheath  radius  to  give  approximate  sheath  radii  for  the 
various  beam  currents  considered.  The  probe  is  found  to  b"  inside  of  the  sheath 
except  for  the  beam  current,  =  1  pA. 

The  potential  of  tne  rear  section  is  shown  to  be  close  to  zero.  It  is  chosen 
equal  to  zero,  which,  witn  the  potential  differences,  determines  the  absolute 
potentials.  An  investigation  of  now  close  V ^  is  to  zero  would  be  interesting.  It  is 
slightly  negative  (positive)  for  ion  (electron)  beams. 

The  total  net  sheath  charge  is  the  negative  of  the  sum  of  the  charges  on  the 
three  rocket  conductors.  Formulae  are  developed  to  give  the  sheath  potential  as  a 
function  of  this  charge,  the  sheath  radius,  and  the  radius  of  the  point  in  the  sheath. 
The  sheath  charge  density  is  assumed  constant.  These  formulae  are  used  to  give 
the  part  of  the  conductor  potentials  that  is  due  only  to  the  sheath.  The  latter  and 
the  conductor  charges  are  determined  by  a  fast  converging  reiteration.  Except 
for  1^  =  1  pA,  the  sheath  potentials  are  relatively  small.  Nevertheless,  use  of 
Jaycor’s  more  realistic  simplified  Lam  sheath  treatment*  instead  of  assuming  con¬ 
stant  charge  density  would  be  interesting  and  not  very  difficult,  to  ascertain  that 
this  more  correct  treatment  would  make  only  a  small  correction. 

The  charge  on  the  rear  section  is  always  positive  for  a  positive  ion  beam, 
negative  for  an  electron  beam,  and  of  smaller  magnitude  than  that  on  the  forward 
section.  The  rear  to  forward  section  voltmeter  current  and  any  probe  emitter 
electron  current  to  the  rear  section  have  a  negligible  effect. 

The  extended  probe  has  a  much  smaller  area  and  is  farther  from  the  forward 
section  than  the  rear  section  is,  so  it  collects  much  less  return  current  for  a  given 
potential.  This  and  the  size  of  the  probe  to  forward  section  voltmeter  resistance 
cause  to  be  between  0  and  V^,  except  for  beam  currents  far  beyond  saturation, 
the  minimum  value  of  1^  for  which  |  vj  =  E^,  the  beam  particle  energy.  Thus, 
in  general,  |<pj  <  1^1 .  This  resistance  is  so  low  that,  for  all  the  positive  ion 
beam  (1^  >  0)  cases  considered,  the  charge  on  the  probe,  q^,  must  be  negative  lor 
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the  probe  and  emitter  current  to  balance  tiie  voltmeter  current.  For  >  0,  much 
larger  values  of  this  resistance  and/or  would  cause  to  be  positive.  For  |  1^1 
beyond  saturation  and  increasing,  tins  voltmeter  current  becomes  more  nearly 
negligible  and  |\^l  approaches  zero.  With  the  voltmeter  used,  as  1 1^1  increases 
from  zero,  at  first  I  V  i  increases  and  is  a  positive  number  appreciably  less 

than  one.  After  saturation  is  approached,  i  V.^l  reaches  a  maximum  then  approaches 
0  (except  for  thermal  energy)  while  0  0^  approaches  1.  Tile  numbers  in  Tables  2 

md  ,1  are  consistent  with  these  conclusions,  but  more  extensive  data  to  test  them 
would  be  desirable,  dunce  the  electron  beam  case  is  way  beyond  saturation,  the 
probe  to  forward  section  voltmeter  resistance  is  effectively  infinite. 

'f  q  =-  0,  the  probeless  potential,  that  is,  the  potential  at  the  position  of  the 
probe  with  (fie  probe  removed,  is  essentially  \  ...  Two  non-equivalent  methods  are 
given  for  calculating  the  approximate  probeless  potential  when  q.^  4  0. 

For  a  positive  ion  beam,  secondaries  reduce  the  sheath  size,  decrease  the 
charge  density  near  the  rocket,  q  ,  and  Iq^l.  f  or  an  electron  beam,  many 
secondaries  may  be  produced  wfiieti  increase  the  electron  density  near  the  forward 
section.  Their  effect  is  to  mcrease  the  electric  field  near  the  rocket,  to  decrease 
tlie  field  m  the  outer  part  the  sheath,  to  increase  qt  and  'q  !,  and,  if  below 
saturation,  to  increase  slightly  the  rocket  potential.  Further  investigation  would 
be  into  re  ting.  The  effect  of  scattering  is  thought  to  be  small  but  should  be 
investigated. 

Ionization  is  negligible  for  positive  ion  beams,  but  should  be  investigated  for 
electron  beams,  for  which  it  lias  a  small  or  moderate  effect. 

The  charging  time  constant  of  the  rocket  (the  forward  section)  due  to  its 
capacitance  is  calculated  for  each  value  of  1^.  It  is  much  larger  than  the  ion  plasma 
period  for  ion  beam  currents  less  than  a  few  50  s  of  p A.  It  dominates  the  time 
behavior  below  ~  10  p  A  ion  beam  current.  The  ion  transit  time  through  the  sheath 
becomes  important  at  higher  beam  currents  and  dominates  the  time  behavior  at  ion 
and  electron  beam  currents  near  or  above  saturation. 

A  number  of  approximations  are  used  in  this  report.  Wore  accurate  treatments 
would  lead  to  appreciable  changes  in  numerical  results,  but  they  would  be  expected 
to  make  little  or  no  qualitative  changes.  A  combination  of  the  capaeitative  coupling 
equations  used  here  and  a  more  complete  probe  theory  may  lead  to  a  prediction  of 
the  potentials  and  charges  on  the  conductors  of  a  rocket  or  satellite  such  as  that 
considered  here.  The  treatment  here  for  three  conductors  may  easily  be  extended 
to  any  number  of  isolated  conductors  or  small  insulators  of  a  rocket  or  satellite. 

From  the  results  of  Lam-Jaycor*  and  this  report,  one  may  induce  the  following 
general  steady  state  properties  of  vehicles  in  the  ionosphere  at  potentials  small 
enough  not  to  cause  breakdown  and  appreciable  ionization  of  background  neutrals. 

For  a  given  beam  current  and  shape  of  vehicle  with  an  all-conducting  surface,  as 
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the  size  increases,  the  sheath  thickness,  vehicle  potential,  and  charge  decrease. 
For  sheaths  that  are  thick  compared  to  vehicle  size,  the  charge  decreases  very 
little.  For  a  given  beam  current  and  vehicle,  as  the  fraction  of  the  surface  that 
is  covered  with  dielectric  is  increased,  the  potential  and  charge  on  the  conducting 
part  increase. 

After  this  report  was  written,  mv  attention  was  drawn  to  a  paper  by  Katz 
8  1 
and  Mandell,  who  made  some  analysis  of  the  rocket  of  Cohen  et  al.  Thev 

concluded  from  their  calculations  that  the  rear  section  is  no  more  than  1  V  negative 

when  the  forward  section  is  charged  negatively  to  several  kV  in  a  plasma  with  a 

fraction  of  an  eV  temperature.  They  indicated  that  this  potential  would  be  positive 

and  of  larger  magnitude  when  the  forward  section  is  positively  charged  over  a  kV, 

but  did  not  consider  the  case  of  its  potential  -  +100  V.  They  also  concluded  that  the 

forward  section  potential  was  nearly  saturated  at  -2  kV  for  1^  =  .174  p  A  (-400  pA). 

Their  explanation  of  the  difference  in  probe  and  rear  section  potentials,  however, 

disagrees  with  that  in  the  present  report.  They  seem  not  to  have  considered  the 

important  effects  of  the  charge  on  the  probe  and  the  probe  to  forward  section 

voltmeter  current. 


8.  Katz,  I.,  and  Mandell,  M.  J.  ( 1982)  Differential  charging  of  high-voltage 

spacecraft:  the  equilibrium  potential  of  insulated  surfaces,  J.  Geophys.  lies. 

87:4533.  -  - 

aw 


23 


References 


1.  Cohen,  H.  A.  ,  Sherman,  C.  ,  and  Mullen,  E.G.  ( 1979)  Spacecraft  charging  due 

to  positive  ion  emission:  an  experimental  studv,  Geophvs.  Res.  Letts.  6:5  15. 

-  v* 

2.  Page,  L.  (1925)  Introduction  to  Theoretical  Physics,  Second  ed.  ,  §120. 

3.  Smvthe,  W.  R.  (  1962)  Charged  right  circular  cylinder,  .1.  Appl.  Phvs.  33:2966. 

-  WVW 

4.  Huber,  W.  8.  private  communication. 

5.  Stannard,  P.  R.  ,  Katz,  I..  Mandell,  M.  .1.  ,  Cassidy,  Parks,  D.  E.  , 

Rotenberg,  M.  ,  Steen,  P.G.  (1980)  Analysis  of  the  charging  of  the 
SCATHA  (P78-2)  satellite.  Report  NASA  CK-165348,  SSS -R -8  1  -47  98, 
pp  22-24,  29.  '  ' . 

6.  Sternglass,  E.J.  (  1954)  Backscattering  of  kilovolt  electrons  from  solids, 

Phys.  Rev.  95:345. 

- - -  MM 

7.  Leadon,  R.  E,  ,  Woods,  A.J.,  Wenaas,  E.  P.  ,  and  Klein,  H.  H.  ( 198 1 )  Analytical 

Investigation  of  Emitting  Probes  in  an  Ionized  Plasma,  A  FG  L -T R  - S  1-0138, 

AD  AWIflff: - - 

8.  Katz,  I.,  and  Mandell,  M.J.  (1982)  Differential  charging  of  high-voltage 

spacecraft:  the  equilibrium  potential  of  insulated  surfaces,  J.  Geophvs.  Res. 
87:4533. 


25 


I 


phboshjo  j*s 


-mot  nus 


Appendix  A 

Auxiliary  Detailed  Calculations 


VI.  FORWARD  TO  RKAR  SKCTION  CHARGF.  RATIO 


In  this  section,  the  ratio  A  of  the  charge  on  the  forward  section  of  the  rocket 
to  that  on  the  rear  section  when  the  two  sections  are  electrically  connected  is  cal¬ 
culated, 

3 

From  Snivthe,  the  surface  charge  density  on  the  side  and  end  respectively  are 
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Figure  Al.  Surface  Charge 
Density  on  Rocket 
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Let  7.  =  by  .  Then  Eq.  (A5)  becomes 
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The  first  integral  is  a  beta  function, 
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These  gamma  functions  are  known  or  may  easily  be  determined  accurately. 

2^  =  0.  41  and  b  =  b^  =  1.475  m  . 
so 

^^2  2 

(1  -  -j=-)  =  0.  521333  .  (A8) 

The  second  integral  in  Eq.  (A6)  was  evaluated  by  expanding  the  binomial,  integrating 
term  by  term,  and  keeping  as  many  terms  as  necessary  for  accuracy.  Since 
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When  z  =  0,  K  =  a,  k  =  1.  777,  and  til  2  '  b  Hu,  then  Q/V  is  within  1.  25%  of  the 
correct  capacitance. 

p  *  is  Q / V  obtained  from  Eq.  (All)  with  a  =  0.  19  m,  k  =  1.777,  b  =  1.27  m, 
K=  1.7  1  m,  |zl  =  1.01  m,  giving 

p130~l  =  °-  2659  111  '  (A15) 

P12  3  obtained  from  Eq.  (All)  with  a  =  0.  19  m,  k  =  1.777,  b  -  1.475  m, 

K=  1.71  m,  |z|  =  0.805  m,  giving 

p17  1  =  0.  264  nf  .  (A  1 6 ) 

Physical  consideration  indicates  that  this  should  be  larger  than  Eq.  (A15).  This 
discrepancy  is  due  to  some  inaccuracy  in  Eq.  (All).  The  latter  may  be  improved 
by  replacing  b4  by  2b4.  This  leads  to  P13q  1  =  0.  260&  nf  and  p^  3  1  =  0.  261^  nf. 
Use  of  these  values  instead  of  Eqs.  (A15)  and  (A16)  would  make  little  change. 


A3.  SHEATH  RADIUS 


In  ttiis  section,  an  approximate  interpolation  formula  for  sheath  radius  is 
developed. 

Leadon  et  al‘  used  a  modification  (simplification)  of  Lam's  method  to  obtain 
the  radius,  rg.  of  the  spherical  sheath  around  an  emitting  probe.  They  equated 
the  thermal  current  density  times  the  sheath  area  to  the  beam  current.  This  gives 

(A  17) 


Fr 

s  V  4  ^Kf  noVt 


I  is  the  return  current  usually  taken  equal  to  I,  ,  q  is  the  charge,  v.  the  modified 
r  /~kT  D  T 

thermal  speed  y-g-  —  —  .  and  m  the  mass  of  the  return  particles. 

The  ion  to  electron  mass  ratio  is  taken  to  be  46,  000.  k  is  Boltzmann's  constant, 

the  plasma  temperature,  T,  is  taken  to  be  550°K,  and  the  number  density,  n  ,  to 

9-3  ° 

be  10'  m  .  Eq.  (A17)  is  valid  when  rg  is  large  compared  to  probe  (rocket)  size. 

When  the  outside  of  the  sheath  is  close  to  the  rocket,  the  sheath  is  assumed  to 

be  a  cylinder  concentric  with  the  rocket.  Similarly,  the  return  current  is  equated 

to  thermal  current  density  times  the  sheath  area 
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where  Rg  is  the  sheath  radius  and  B  is  half  the  sheath  length.  The  latter  is 
assumed  to  be 


B  =  (b/a)R  . 

s 


(A  19) 


So  the  sheath  radius  is  assumed  to  be 
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For  small  beam  currents,  f  =  ;  for  large  beam  currents,  f  =  j. 

The  interpolation  factor  is  taken  to  be 

(x  +  I  )a 
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Somewhat  arbitrarily,  x  is  chosen  so  that  f  =  f ^  =  f^fO.  5)  s  — )  when 

R„  =  b.  Eq.  (A20)  is  solved  for  1=1.  when  R  =  b  and  f  =  f.  .  Then  Eq.  (A21)  is 
s  r  rb  s  b 

solved  for  x  when  f  =  f.  and  I  =  I  .  .  This  leads  to 
b  r  rb 
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for  an  ion  beam,  and 
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for  an  electron  beam,  where  R  is  in  meters  and  I  in  a  A.  The  values  of  R  in 

s  r  s 

Table  3  are  calculated  from  Eqs.  (A22)  and  (A23).  Values  are  also  calculated 

the  same  way  but  with  f.  =  f.  (0.  7)  =  0.  3  X  +  0.  7  X  - . 

b  b  £•  2b  +  a 


A4.  THE  POTENTIAL  DUE  TO  THE  SHEATH 

For  simplicity,  the  sheath  charge  density  p,  is  assumed  constant  in  this  first 
treatment  of  the  sheath. 

First,  suppose  a  spherical  conductor  of  radius  a  and  charge  -q  surrounded  by 
a  spherical  sheath  of  radius  rg  and  total  charge  qg  (Figure  A2). 
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Figure  A2.  Spherical 
Conductor  and  Sheath 


From  Gauss's  Law,  for  a  5  r  <  rg. 

4  n  r2E  =  4  n  [-qg  +  ~  ?r(r3  -  a3)p]g  , 


(A25) 


where  E  is  the  (radial)  electric  field  at  a  radial  distance  r.  The  potential,  V,  at 
r  minus  the  potential  there  due  to  the  conductor  is  the  potential  due  to  the  sheath. 
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Secondly,  the  potential  is  obtained  in  the  median  plane,  z  =  0,  due  to  a  cylindrical 
sheath  of  radius  Rg  and  half  length  B  -  (b/a)Rg  of  charge  qg  surrounding  a  cylindrical 
conductor  of  radius  a  and  half  length  b  (Figure  A3).  First,  the  case  of  Rs  R  is 
treated.  Eq.  (All)  is  used  with,  for  simplicity,  z  =  0.  The  potential  at  R  due  to 
an  infinitesimal  shell  at  R1  £  R  and  z'  £  z  is 
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Figure  A3.  Cylindrical 
Conductor  and  Sheath 


The  potential  at  R  due  to  a  similar  shell  at  R'  2  R  and  21  -  z  is  assumed  independent 
of  R,  so  is 
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To  make  the  exponent  constant  and  thus  the  potential  integrable.  z'  is  taken  =  (b/a)R'. 
The  charge  density  is  qg  divided  by  the  sheath  volume. 
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Thus,  the  potential  due  to  the  sheath  at  a  point  R  inside  the  sheath  is 
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This  reduces  to  Eq.  (A26)  for  b  =  a  and  k  =  1.  For  the  case  R  -  R  ,  using  Eq.  (A27), 
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As  R  -*  *>,  V  -*  — as  it  must, 
sc  R 

For  increasing  Rg  >>  b,  the  sheath  surface  becomes  spherical.  So,  the  poten¬ 
tial  due  to  such  a  sheath  at  point  R  in  the  median  plane  is  calculated.  It  is  the 
potential  due  to  the  sphere  of  uniform  charge  density  p  minus  that  due  to  the 
cylinder  a,  b  of  density  p.  First,  this  potential  for  R  s  Rg  is  calculated. 


i  Ra3pg  f3  T*  R'3pg 

'  ss(R>  =  1— W +  J  - 2 dR' 

88  Ks  R  R' 
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Figure  A4.  Cylindrical 
Conductor  and  Spherical 
Sheath 


is  given  by  Eq.  (A32)  in  the  limit  of  a  -»  0  at  constant  b/a,  then  Rg  is  replaced 
by  a  and  qg  by  ira^2bp. 

( j-JtRg  -  7ta^2b)  p  *  qg  .  (A34) 


So,  Eq.  (A33)  becomes 


V*  JR)  = 

S3 


qsg 

3 — n: 

•773  b 


(A35) 


X 


rkab 


b(kF  +  1} 


(A36 ) 

If  k  =  1  and  b  =  a,  this  reduces  to  Eq.  (A26)  with  a3  replaced  by  (3/2)a3.  If  k  =  3/2 
andb=(2/3)a,  Eq.  (A35)  reduces  to  Eq.  (A26).  Secondly,  the  potential  due  to  the 
sheath  at  R  2  Rg  is 


,(R)  = 


Rg 


P  g 


-  V  = 
c 


V 


Tir 

?a  b 


R  3 
s 


-  X 


(A37) 
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As  i-  -*  ■*.  V  -»  — tv  as  it  must, 
ss  K 

The  sheath  potential  is  taken  as 


V  (K>  =  F V  +  (1  -  F)V 
s  ss  sc 


(A38) 


both  for  H  -  B  and  H  -  R  .  F  -»  1  as  R  -»  -r  and  F  -*  0  as  R 
s  s  s  s 

interpolation  factor  is  taken  to  be 


R  -  a 
s 

Tr~T7 


So,  the 


0.3C0 


Somewhat  arbitrarily,  y  is  chosen  so  that  F'  =  F,  =0.5  when  R  =  b.  Then, 

■  b  s 

v  =  b  -  2a.  The  same  calculation  is  also  made  with  f.  =  f.  (0.  7)  and  F,  =  0.  7,  for 

b  b  b 

which  y  =  (3b  -  10a)/ 7. 


A5.  ION  TRAVERSE  TIME  THROUGH  THE  SHEATH 


For  simplicity,  the  traverse  time  is  taken  to  be  the  time  calculated  for  an  ion 
to  move  radially  through  a  sheath  around  a  charging  spherical  conductor,  starting 
when  the  beam  is  turned  on  at  t  =  0.  The  beam  current  rise  time  is  assumed  to  be 
much  shorter  than  this  traverse  time.  Also,  the  assumption  is  made  that  the  field 
grows  exponentially  with  time  and  the  radial  dependence  of  field  and  potential  art- 
independent  of  time.  The  radial  acceleration  of  an  ion  of  charge  e  and  mass  M  is 


e 

M 


(A40) 


The  values  of  e  and  M  taken  are  the  magnitude  of  the  electron  charge  and  the  mass 
respectively  of  an  NO+  ion.  From  Eqs.  (A24),  (A25),  (A40)  and  the  time  assumption. 


eV 

M ( r  3  -  a3) 
s 


(A41) 


The  value  of  rg  is  taken  from  row  3  of  Table  3.  Integrating  E  from  a  to  rg  at 
t  =  oo  leads  to 


V(a> 


r^  -  (r  -  a)a/2 
s  s 

— 2 - r 

r  +  r  a  +  a 
s  s 


(A42) 
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The  value  of  a  is  chosen  so  that  Via)  equals  the  value  in  row  4  of  Table  3.  This 
results  in  0.  47  <  a  <  0.  51  m. 

For  an  ion  beam,  at  t  =  0;  r  =  r  and  dr/dt  =  -v. .  From  Section  A3 

=  156  m/s.  The  mass  is  taken  to  be  that  of  NO  .  The  traverse  time,  t^,  is  the 

time  at  which  r  =  a.  For  the  electron  beam,  at  t  =  0:  r  =  a  and  dr/dt  is  chosen  =  0. 

The  traverse  time  is  the  time  at  which  r  =  r  .  The  value  of  t.  is  determined  either 

s  t 

by  trial  and  error,  solving  Eq.  (A41)  numerically,  or  by  solving  numerically 


r  -.Q  r  3  t 

*  3  r  l-e'~ 

j  [  j 


The  results  for  r  1  =  7  and  for  7  '  =  >  T  are  shown  in  rows  19  and  20  respectively 

of  Table  3. 


